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Abstract 

We apply the effective field theory approach to the coupled metric-inflaton system, in order to 
investigate the impact of higher dimension operators on the spectrum of scalar and tensor pertur- 
bations in the short-wavelength regime. In both cases, effective corrections at tree-level become 
important when the Hubble parameter is of the order of the Planck mass, or when the physical wave 
number of a cosmological perturbation mode approaches the square of the Planck mass divided by 
the Hubble constant. Thus, the cut-off length below which conventional cosmological perturbation 
theory does not apply is likely to be much smaller than the Planck length. This has implications 
for the observability of "trans-Planckian" effects in the spectrum of primordial perturbations. 
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I. INTRODUCTION 



One of the main successes of inflation [1] is the explanation of the origin of structure [2]. 
During slow- roll, the Hubble radius remains nearly constant, while cosmological modes are 
constantly pushed out of the horizon. Thus, local processes determine the amplitude and 
properties of perturbations at sub-horizon scales, which are transfered to cosmologically 
large distances by the accelerated expansion. In that sense, the sky is the screen upon which 
inflation has projected the physics of the microscopic universe. 

The primordial perturbations seeded during inflation arise from quantum-mechanical fluc- 
tuations of the inflaton around its homogeneous value. Hence, their properties directly de- 
pend on the quantum state of the inflaton perturbations. Conventionally, this is taken to be a 
state devoid of quanta in the asymptotic past, raising the crucial question of whether we can 
trust cosmological perturbation theory — and its quantum nature — at such early times [3]. 

According to our present understanding, quantum field theories and general relativity 
are merely low energy descriptions of a more fundamental theory of quantum gravity. In 
the case of inflation, the leading terms in the corresponding effective Lagrangian are the 
Einstein-Hilbert term plus the inflaton kinetic term and potential. In an effective field 
theory treatment, these terms are accompanied by all other possible operators compatible 
with the symmetries of the theory, namely, general covariance and any other symmetry of 
the inflaton sector. Higher dimensional operators are suppressed by powers of an energy 
scale, which we will assume to be of the order of the reduced Planck mass M p = (8-7rG) _1 / 2 , 
and they are therefore expected to be negligible at sufficiently small momenta, or sufficiently 
long wavelengths. The goal of this paper is to determine the three-momentum scale A at 
which these operators become important. The terms that yield the leading (momentum- 
independent) corrections to the primordial spectrum have been recently discussed in [4]. 

It is crucial to realize that A does not need to equal the Planck scale. As an extreme 
example, let us consider perturbations around a background given by a constant scalar field 
in Minkowski spacetime, (p = ip + 5ip and = rj^ + Sg^. In this case, the effective action 
for the perturbations is Lorentz invariant, because the background is Lorentz invariant and 
we assume that the original effective action S[<p, g^ v ] is generally covariant (we are implicitly 
assuming that a Lorentz- invariant gauge on 8g^ u has been imposed.) Therefore, the regime 
in which higher order operators become important will be expressed by Lorentz invariant 
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relations such as, for example, k^ « Mp. In other words, there is not any bound on 
the three-momentum |k|, simply because a relation like |k| « M p is not Lorentz invariant. 
One would like to extend this argument to an expanding universe. On short time-scales and 
short distances, spacetime can be regarded as flat and the scalar field driving inflation can be 
regarded as constant. Hence, it seems that cosmological perturbation theory should be valid 
all the way into the regime |k| — > oo, because the relevant four-momenta of cosmological 
perturbations are light-like, k^k^ = 0. As we shall see though, the evolution of the inflaton 
leads to small but finite violations of the Lorentz symmetry even in the short-wavelength 
limit, which are imprinted on the power spectrum as k-dependent corrections. 

The phenomenological imprints of trans-Planckian physics on the primordial spectrum of 
perturbations, and the implications of a finite cut-off A on the spatial momentum of cosmo- 
logical modes have been extensively studied [5-30]. These articles mostly study corrections 
to the power spectrum in the long- wavelength limit k/a = k p h <C H, at late times, which 
is the regime directly accessible by experimental probes. In this article we focus instead on 
the short-wavelength regime k p h ^> H, at early times, since we are interested in determining 
how far into the ultraviolet cosmological perturbation theory applies. At short wavelengths, 
the power spectrum can be cast as a derivative expansion of the form 

<^*(k)Mk)> = ^ (l + « 2 || + + • • • ) , (1.1) 

with coefficients a, that depend on slow-roll parameters and the dimensionless ratio H/M p . 
The analytic corrections to the leading result l/2k arise from tree-level diagrams with ver- 
tices from higher-dimensional operators. We only consider tree-level diagrams here, since 
we expect loop diagrams to simply introduce a logarithmic dependence of the dimensionless 
coefficients ccj on scale. Cosmological perturbation theory fails when the expansion in powers 
of k p h breaks down, namely, when all the terms become of the same order, 

k ph ^M p .f^ = A. (1.2) 

As we shall show, the ratios a 2 n/«2n+2 are all quite large and of the same order, so the 
effective cut-off A significantly differs from M p . In a slightly different context, a similar 
analysis has been applied to the bispectrum in [31]. 

The structure of this paper is as follows. In the next section we describe the relevant 
background to our problem, setting up a description of perturbation theory in cosmology and 
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describing the basics of tensor and scalar fluctuations. In section III we compute the squared 
amplitude of tensor perturbations and we derive the results mentioned above. In section IV 
we apply a similar analysis to the case of scalar perturbations, and obtain similar results. 
We conclude and discuss possible implications of our results in section V. Throughout this 
paper we work in units such that h = c=l, and with metric signature ( — h ++)• 



II. SETTING THE SCENE COSMOLOGICAL PERTURBATION THEORY 
A. The Inflating Background 

Our starting point is a standard single-field inflation model. At sufficiently late times, 
the inflaton and gravity must be described by a low-energy effective action, whose leading 
terms are dictated by general covariance and the field content, 



S = d 4 x^g 



Ml 1 

-±R--dp<p&*<p-V{<p) 



(2.1) 



In an effective field theory context, the action should also contain additional terms sup- 
pressed by powers of a dimensionful scale, which we assume to be of the order of the Planck 
mass M p . Our goal is to determine the point beyond which such higher-dimensional op- 
erators produce corrections to the two-point function of cosmological perturbations that 
cannot be neglected. Our considerations can be readily generalized to cases in which the 
suppression scale of the higher- dimensional operators is not the Planck mass, but any other 
scale. 

If the potential V(ip) is sufficiently flat, at least in a certain region in field space, there 
exist inflationary solutions, along which a homogeneous scalar field ip(rj,x) = (fo(rj) slowly 
rolls down the potential and spacetime is spatially homogeneous, isotropic and flat, 

9$ = * 2 (V)V»», (2-2) 

where i]^ is the Minkowski metric and r\ denotes conformal time. A model-independent 
measure of the slowness of the inflation is given by the slow-roll parameter 

where H = a' /a 2 is the Hubble parameter and a prime denotes a derivative with respect 
to conformal time. During slow- roll, e is nearly constant, and to lowest order in slow- roll 



parameters its time derivative can be neglected. Throughout this article we work to leading 
non-vanishing order in the slow-roll expansion. 

B. Cosmological Perturbations 

Let us now consider cosmological perturbations around the homogeneous and isotropic 
background described above. Writing tp = (p + dp and g^ = g$(ri) + Sg^r), x), and 
substituting into equation (2.1), we can expand the action So up to the desired order in the 
fluctuations 5p and Sg^ u , 

Solp^g^} = 5 S + 5 1 S + 5 2 S ^ . (2.4) 

The lowest order term SqSq does not contain any fluctuations and describes the inflating 
background; the linear term 5\So vanishes because it corresponds to the first variation of 
the action along the background solution, and the quadratic part of the action 82 So describes 
the free dynamics of the perturbations. The latter is what we need in order to calculate the 
primordial spectrum of fluctuations. To quadratic order, tensor and scalar perturbations 
are decoupled, so we may study them separately. 

1. Tensor Perturbations 

Tensor perturbations are described by a transverse and traceless tensor hy, 

ds 2 = a 2 (r]) [-drj 2 + (<% + h i j)dx i dx i ] . (2.5) 

The tensor hij itself can be decomposed in plane waves of two different polarizations, which 
again decouple at quadratic order. We shall hence focus on just one of them, 

w k 

where the hk(r]) are the corresponding mode functions, and e^k) denotes the normalized 
graviton polarization tensor, e^-e^j = 1 (we raise and lower spatial indices with the Kronecker 
delta.) Note that for later convenience we work in a toroidal universe of volume V = L 3 ; 
hence, the spatial wave numbers have components ki = rii- (2-k/L), where the rij are arbitrary 
integers. 
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Substituting the expansion (2.6) into the action (2.1), and using the background equations 
of motion, we may then express the free action (^Sq as 



\ / *' E 

J 1, 



v' k v'_ k - [ k 2 - 



V k V-k 



(2.7) 



where the scalar variable v k is defined as 



(2.8) 



Thus, in terms of v k the action for tensor perturbations takes the form of an harmonic 
oscillator with time- dependent frequency. 



2. Scalar Perturbations 



In spatially flat gauge, the perturbed metric reads 

ds 2 = a 2 (r]) [-(1 + 2(p)dr] 2 + 2d i Bdx i dr ] + c^dxW] . (2.9) 

On first inspection there appear to be three independent scalar variables: <p, B and the 
inflaton perturbation 5ip. However, Einstein equations impose constraints on both and B. 
Solving the corresponding Fourier transformed equations to leading order in the slow-roll 
expansion, one finds (see e.g. [32]) 

Consequently, there is only one physical scalar degree of freedom, and scalar perturbations 
can be described by just one variable. A convenient choice that is particularly useful for 
quantizing scalar perturbations is the Mukhanov variable [33], which in spatially flat gauge 
takes the simple form 

v k = a5ip k ■ (2-H) 

Using relations (2.10) and (2.11), we may express S2S0 in terms of v k only. For constant 
e, that is, to leading order in the slow-roll expansion, the resulting action is also given by 
equation (2.7). This agreement greatly simplifies the analysis, because it allows us to use 
the same set of propagators to describe both scalar and tensor fluctuations. 
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To leading order in the slow-roll expansion, the mode functions of both scalar and tensor 
perturbations hence satisfy the same equation of motion during inflation. Varying the action 
(2.7) with respect to v_ k we obtain 

v k = 0, (2.12) 

which has a unique solution for appropriate initial conditions. The conventional choice is 
the Bunch-Davies or adiabatic vacuum, whose mode functions obey 

l + o(±)]. (2.13) 

Because we are only interested in the sub-horizon limit, this is all we need to know about 
the mode functions. In particular, because the behavior of the mode functions in the short- 
wavelength limit does not depend on the details of inflation, our results are also insensitive 
to the particular form of the inflaton potential. 



k 2 - 



w) — ► 



2k 



C. Quantum Fluctuations and the in-in Formalism 

In order to study the properties of cosmological modes in the short-wavelength regime, 
we concentrate on the two-point function of the field v, 

(v*(Tj,k)v(ri,k)) = (0,in\v*(r],'k)v(ri,k)\0,in) , (2.14) 

where \0,in) is the quantum state of the perturbations, which we assume to be the Bunch- 
Davies vacuum. The two-point function characterizes the mean square amplitude of cos- 
mological perturbation modes, and differs from the power spectrum just by a normaliza- 
tion factor. Note that in an infinite universe, the two-point function is proportional to a 
momentum-conserving delta function, which in a spatially compact universe is replaced by 
a Kronecker delta. 

In the in-in formalism (see [34] for a clear and detailed exposition) the two-point function 
can be expressed path integral, 

(u*(77,k)v(77,k))= (2.15) 
= J Vv + Vv„ vl(r), k)v_(ri, k) exp (iS irce [v + , v_)) exp (iS int [v + ]) exp {-iS int [v_]) , 

where Sf ree is quadratic in the fields, and S- m t contains not just the remaining cubic and higher 
order terms in the action, but also any other quadratic terms we may decide to regard as 
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perturbations. Note that there are two copies of the integration fields v_ and v + , because we 
are calculating expectation values, rather than in-out matrix elements. This path integral 
expression is very useful to perturbatively expand the expectation value in powers of any 
interaction. In particular, each contribution can be represented by a Feynman diagram, 
with vertices drawn from the terms in Si nt and propagators determined by the free action 
Sfrec- In our case, the latter are given by 

e -ik\r)-r)'\ 

Vv + Vv_ v* + (r], k)v + (?/, k) exp(zS , frcc ) « — — — , (2.16a) 



T) T}' 
• • 

rj rj' 



■ • 



Vv + Vv_ v*_(v, k)v_(r)', k) exp(^ frce ) « ^l^i (2.16b) 

2k 

e ik{ri~ri') 

Vv+Vv_ v* + (r], k)v_(r)', k) exp(iS hcc ) « — — — , (2.16c) 



rj rj 1 

which we quote here just in the sub-horizon limit. Note that to first order in Si nt there 
are two vertices, one that contains powers of v + and one that contains powers of V-\ the 
associated coefficients just differ by an overall sign. 1 

As a simple example, let us calculate the value of the two-point function in the short- 
wavelength limit to zeroth order in the interactions. Using the definition (2.15) and equation 
(2.16c), we find 

(v*( V , k)v( V , k)> = . • « (\k V \ > 1), (2.17) 

TJ TJ AK 

which is the well-known and standard short-wavelength limit result. In this regime, the 
two-point function is hence the inverse of the dispersion relation, since the latter determines 
the appropriate boundary conditions for the mode functions, as in equation (2.13). 

In the next two sections we use the path integral (2.15) to calculate the corrections to 
the two-point function coming from higher-order operators in the action. These can be 
interpreted as corrections to the dispersion relation, even though in the presence of such 
terms the mode equations generally contain higher order time derivatives. In any case, 
a significant disagreement between the calculated two-point function and the lowest order 



1 The quadratic action Sf r ee enforces w+(fc) = u_(fc) at time rj. Hence, we could replace u* (rj, k)v_ (77, k) 
by v* (ji, k)i> + (?7, k) or v* (77, k)w_ (77, k) inside the path integral (2.15). Our choice removes the apparently 
ill-defined corrections we otherwise obtain when higher-order time derivatives act on the time-ordered 
products in equations (2.16a) and (2.16b). These ill-defined corrections can also be eliminated by field- 
redefinitions, a procedure that leads to the same corrections we find using our choice of field insertions. 
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result (2.17) points to the lack of self-consistency of our quantization procedure, and signals 
the breakdown of cosmological perturbation theory. 

III. THE LIMITS OF PERTURBATION THEORY: TENSORS 

The lowest order action (2.1) contains the leading terms that describe the dynamics 
of the inflaton and its perturbations. However, as we have noted, in an effective field 
theory approach the action generically contains all possible terms compatible with general 
covariance and any other symmetry of the theory. Here, for simplicity, we assume invariance 
under parity, an approximate shift symmetry of the inflaton, and a discrete Z 2 symmetry 
(f — > —tp. Thus, all possible effective corrections to the action (2.1) can be built from the 
metric g^ u , the Riemann tensor R^xp, the covariant derivative V M and an even number 
of scalar fields tp. In what follows, we consider these additional terms and compute the 
corrections they induce on the two-point function of tensor perturbations in the short- 
wavelength limit. This allows us to determine the regime in which additional terms in 
the action cannot be neglected, and hence, the range over which cosmological perturbation 
theory is applicable. The reader not interested in technical details may skip directly to 
section III C, where we collect and summarize our results. 

A. Dimension Four Operators 

We begin our analysis by considering all dimension four operators, which appear in the 
action multiplied by dimensionless coefficients. On dimensional grounds, we expect these 
to yield corrections to the two-point function that are suppressed by only two powers 2 of 
M p . These operators will also help us to illustrate our formalism and discuss some of the 
important issues related to our calculation. 

Any generally covariant dimension four effective correction must be of the form 




(3.1) 



where C 2 



is the square of the Weyl tensor, 



— D , T3liv\p 



2R lil/ R ,MU + g-R 2 ; 



(3.2) 



2 Dimension six operators quadratic in ip also contribute at this order; we consider them later. 
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and the dimensionless couplings a and (3 are assumed to be of order one. Note that we 
have ignored total derivatives like the Gauss-Bonnet term, since they do not lead to any 
corrections in perturbation theory. The Levi-Civita tensor cannot appear in the action 
because we assume invariance under parity. 

We start by substituting the perturbed metric (2.5) into equation (3.1) and expanding 
up to second order in hij. Using the modified background equations and equation (2.8) to 
express the tensor perturbations in terms of the variable v, we obtain in the sub-horizon 
limit 

= ^E J 'dn\~v k [vU + tfv- k }-^[vl + k\ t ]v- k y (3.3) 

- ^E / i/jlK^i-^^j.jiK^-M^)}. 

From these expressions, it is easy to derive the rules for the vertices 
a a 

= X (3.4a) 

— *— * W iyv'{(^ + k 2 )-^^aH}^{^ v , + e)-2aH-d v/ } 

= x , (3.4b) 

where the arrows indicate the propagator on which the derivatives act (because the vertex 
is quadratic, two propagators meet at the vertex.) 

We are now ready to consider the correction due to the square of the Ricci scalar. The 
first order correction to the two-point function is given by the sum of the following two 
graphs, 

a ia P , . (i5(ri-ri')Qa") a 12H 2 

-X • « — — / dr] ' 



r] r] AftJ-n ' { 2k a 3 j 2k M£ 

a ( a V 
. _ x . = . H . 3.5 

r] V \V vJ 

where we have used the fact that a" /a 3 ~ 2H 2 to lowest order in slow-roll. Notice that the 
operator {^d 2 , + k 2 ^j acting on the time-ordered propagators (2.16a) or (2.16b) produces a 
delta function, since both are Green's functions. On the other hand, when the same operator 
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acts on the propagator (2.16c) we get zero, because the latter is a regular solution of the 
free equation of motion (2.12) in the sub-horizon limit. This remark will turn out to be very 
useful when studying higher dimension operators. 

We can now consider the correction due to the square of the Weyl tensor. In this case, 
the first order contribution is given by the sum of the following two graphs 

I 2i/3 ? J , r_, , u2 o lkM ^\ ^ (3 UiHk ph , 2H 2 



— oo 

■--x • = • x • . (3.6) 



V V w V , 

Note that the imaginary parts cancel once we sum the two graphs. This result is quite 
general and ensures that only corrections with even powers of k p h appear. 

In conclusion, we have found that the leading corrections due to dimension four operators 
result in a two-point function which in the short- wavelength limit has the form 



a 



P \ 1 



• • + • X • + • X • + c.c. « — 

7] 7] \7] 7] 7] 7] J 2k 



H 



p 



(3.7) 

Thus, when H becomes of order M p , these corrections become as important as the leading 
result, and standard cosmological perturbation theory ceases to be applicable, as the reader 
may have expected. 



B. Higher Dimension Operators 

We would now like to consider a generic operator of dimension 2d + 4, suppressed by a 
factor of order 1/M 2d . However, it turns out that considering directly corrections to the 
action (2.7) for the perturbations is a much more efficient approach than starting from 
generally covariant effective terms added to the Lagrangian (2.1), particularly if we are 
interested in identifying the dominant corrections in the sub-horizon limit. Hence, we shall 
focus directly on modifications to the action for the perturbations. A related approach has 
been described in [35]. 

Dimensional analysis implies that any operator of dimension 2d + 4, quadratic in the 
dimensionless tensor perturbations and proportional to 2f powers of the inflaton field (p 
must contain 2d — 2f + 4 derivatives <9 M acting on h^, fo{v) an d a {v)- The derivatives can 
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be distributed and contracted using the Minkowski metric in many different ways, 3 but each 
of these terms can be schematically represented as 

where we have used equation (2.8) and 9* [a, <£>o] is just a symbol that represents any com- 
bination of i derivatives acting on a's and <^o's- One such term would have 2n + m + p 
derivatives acting on one or more factors of a or <p , 2q + m + r derivatives acting on one 
field v and 2s + p + r acting on the other v. In particular, 2n of the derivatives acting on 
the scale factor or the background field are contracted among themselves while m and p of 
them are contracted with derivatives acting on, respectively, the first and second field v. 
The derivatives acting on the fields v are organized in a similar way. 

Let us illustrate this notation by considering a term with p = s = f = 0, m = n = q = 1 
and r = 2. Dimensional analysis implies that d — 3, and thus the explicit form of such a 
term would be 

M p - 8 «9 2+1+0 [a] d 2+l+2 v <9o+o+2 v = M; 8 d^ d v [a] d x d x d v d a d p v d a d? v , (3.9) 

where d^d^d v [a] denotes all possible ways to construct a term with three derivatives of the 
scale factor, with the given tensor structure. 

The first step to estimate the leading correction due to a term of the form (3.8) is realizing 
that this can always be expressed as a linear combination of terms of the form 

M -2d-2 d 2j+l d 2m+l v q^ v ^ ( 3 _ 1Q) 

plus, possibly, a term with no derivatives acting on v, which in any case gives a contribution 
that is always subdominant in the sub-horizon limit. For a proof that this decomposition is 
always possible, we refer the reader to Appendix A. In what follows, we therefore restrict 
ourselves to terms of the form (3.10). 

3 The reader may think that derivatives could be contracted not only among each other with the Minkowski 
metric, but also by using the additional tensor structure provided by the metric perturbations Sg^/a 2 = 
hijVuiVvj- However, it turns out that [Sg^/a 2 ^ d u a — hijij^d^a = and, since hij is transverse, 

d» (Sg^/a 2 ) = /,,., ...T/,, ; = . 

Thus, we get a non-vanishing contribution only when we contract derivatives with the Minkowski metric 
while the factors rj^iTj^j are contracted among each other yielding an irrelevant overall factor. 
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Dimensional analysis requires that the indexes j, I and m in equation (3.10) obey 



j + l + m = d- f + 1 . (3.11) 



Furthermore, since d n (po/dr] n oc V2e M p a n H n and d n a/dr] n oc a n+1 H n , each field ipo yields 
a factor of M p , while each derivative acting on it or on the scale factor results in a factor 
of H to leading order in slow-roll. Finally, the I partial derivatives <9 M acting on v that are 
contracted with derivatives acting on a or tpo can be turned into derivatives with respect to 
7] only. Thus, (3.10) can be re- written as 

M -2d+2f-2f^ H 2j+l D m d l v ^ y ^ ( 3 _ 12 ) 

where we have defined □ = <9 M <9 M ; the corresponding correction to the two-point function is 
schematically given by 

"* ' = ; J- V+ 2 P drj'f{a)H^ 1 x (3.13) 



x . . ( n m *d l v , □ + □ ~d l n , n m ) . . 

rj v' v' V 

plus the complex conjugate of this graph. Because (2.16c) satisfies the free equation of 
motion, this correction is non-vanishing only when the index m is equal to zero, and in this 
case we obtain 

* 1 = / " I'M = (3-14) 



17 77 M™- 2 > +2 J ' "2k 2k Ml i ~ v+2 

-OC 



The leading correction in the short-wavelength limit is the one with the maximum num- 
ber of powers of k. According to equation (3.11), this maximum number simply equals 
d — f + 1 = Imax, and it corresponds to the case in which j = m = 0. Thus, if d — f is 
odd, Imax is even and the leading correction is simply given by 

^.( kW k)> cc JL (iLp' (h±f' +1 (d -f odd), (3,5) 

since each factor of k/M p must be accompanied by a factor of a to render the spatial 
momentum physical. On the other hand, if d — f is even, l max as defined above is odd and 
the term with the highest number of powers of k is purely imaginary. As we have seen in the 
previous section, such a term disappears when we add the contribution from the complex 
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d-f 


Leading correction 


d-f 


Leading correction 





H 2 /Ml 


4 




1 


H 2 klJM* 


5 




2 




6 




3 




7 


H*klJMl« 



TABLE I: Leading corrections to the gravitational wave two-point functions in the short- wavelength 
limit. 

conjugate graph. Therefore, the leading correction corresponds to the largest even value of 
/, which turns out to be l max = d — f, and is therefore given by 



6{v*(k)v(k)) oc 



2k \M 



(d — f even) . 



(3.16) 



Equations (3.15) and (3.16) represent the main results of this section: they express the 
leading corrections to the two-point function (in the sub-horizon limit) associated with a 
generic operator of dimension 2d + 4 containing 2f powers of the inflaton field. Since we 
have assumed an approximate shift invariance of the inflaton, the total number of derivatives, 
2d — 2f + 2, must be greater or equal than the number of fields 2f, which in turn implies 
that d > f. Thus, we can label all the possible corrections with the non-negative index 
d — f. Their magnitude is given in Table I for the first eight values of d — f. Note that 
corrections with d—f — arise from the operators identified by Weinberg in [4]. The leading 
momentum-dependent corrections are given by operators with d — f = 1. 

So far, we have calculated the largest possible corrections to the two-point function in 
the sub-horizon limit given a certain value of d — f. However, the reader might still wonder 
whether such terms can be actually obtained from a covariant action. Employing the same 
technique we used to study the impact of the lowest order terms, it is indeed possible to 
show — after some rather lengthy calculations — that the following family of covariant terms 
generates this kind of contributions, 



d - f = : R^R^ 
d-f = l: (V a R^) V Q iV 
d-f = 2: (V Q V /3 i?^) VaVpR, 



(3.17) 



fJ,V 
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It can be also verified that the d — f = 1 term yields a correction to the two-point function 
proportional to the slow-roll parameter e, and given the structure of this family of operators, 
we anticipate the remaining terms to share the same slow-roll suppression. 



C. The Breakdown of Perturbation Theory 



The corrections to the two-point function are functions of two dimensionful parameters, 
H and k p h- For our purposes, it is convenient to organize these corrections in powers of 
k p h- Thus, following Table I, and reintroducing the subleading terms that we previously 



neglected, we find that the two-point function is 



(v*(k)v(k)) 



2k 



H 2 

±v± p 



H 2 



H 4 



H 4 



H 6 



k 4 
K v h 



k 2 

Jl 2 
p 



+ 



+ 



(3.18) 



The coefficient «2o is of order one, while all the a nn with n > 2 are of order e, as the family 
of covariant terms (3.18) suggests. At the end of Section IV we provide further evidence 
supporting this claim. 

In order for Equation (3.18) to be a valid perturbative expansion, every correction term 
must be much smaller than one. Because o^o is of order one, this implies the condition 



H 



(3.19) 



which must hold for all values of k p h- Equation (3.18) then shows that if condition (3.19) is 
satisfied, the corrections to the two-point function remain small even for k p h ~ M p . In fact, 
to leading order in H/M p we can rewrite equation (3.18) as 



(v*(k)v(k)) « - 



k 2 k^ 
1 + «22^Y + «44^ + 



(3.20) 



where we have introduced the effective cut-off 

A 



H 



(3.21) 



Equations (3.20) and (3.21) are the main result of this article. For k p h A, all the 
corrections are strongly suppressed and can thus be neglected. However, at k p h ~ A, all the 



15 



corrections become of order e, the asymptotic series breaks down, and the effective theory 
ceases to be valid. 

To conclude this section, let us briefly comment on the effects of terms that break the 
shift symmetry. Because the only difference is that these terms contain undifferentiated 
scalars, any such correction can be cast as a generally covariant term that respects the 
symmetry, multiplied by a power of the dimensionless ratio ip/M p . Hence, these terms 
introduce corrections to the two-point function of the form we have already discussed, but 
with coefficients that can now depend on arbitrary powers of the background field (fo, 



Therefore, in the absence of any mechanism or symmetry that keeps the coefficients 
a\j ,a^\ . . . small (e.g. an approximate shift symmetry), such an expansion looses its va- 
lidity for ip Q > M p , regardless of the value of k p h- If, on the other hand, equation (3.22) 



assume that the shift-symmetry is exact, and perturbations theory breaks down again at 
k ph « A. 

D. Loop Diagrams 

Our analysis so far has concentrated on tree-level diagrams, though loop corrections could 
also invalidate the short-wavelength expansion. The contribution of loops is obscured by the 
appearance of divergent momentum integrals. At any order in the derivative expansion it 
is still possible to cancel these divergences by renormalizing a finite number of parameters, 
provided that all terms consistent with the symmetries of the theory are included in the 
Lagrangian [36]. In practice, this cancellation is due to the presence of appropriate counter- 
terms in the Lagrangian. For this reason, divergent integrals in loop diagrams are rather 
harmless. They yield corrections of the same structure as tree-level diagrams, modulo a 
(mild) logarithmic running of their values with scale [37]. Hence, we do not expect this type 
of contributions to drastically change our conclusions, though we should emphasize that this 
is just an expectation. 




(3.22) 



is a sensible expansion, and off is much greater than a^\a^\ then we can effectively 
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IV. THE LIMITS OF PERTURBATION THEORY: SCALARS 



We now turn our attention to corrections to the two-point function of scalar perturbations. 
Despite some complications that are particular to the this sector, the method developed in 
the previous section can be easily extended to scalars. 

To this end, let us consider the action S = So + ASi, where So is the action (2.1) describing 
a scalar field minimally coupled to Einstein gravity, while Si is a generic generally covariant 
correction suppressed by a coupling A ~ 1/M 2d . As we pointed out in the previous section, 
Si generically involves contractions of the Riemann tensor R^xp and the covariant derivative 
V M as well as the scalar field ip. In order to compute the resulting first order contribution to 
the two-point function for v = a Sip, we insert the perturbed metric (2.9) and the perturbed 
inflaton field into the action S and expand up to second order in Sip, 4> and B. We then 
express the metric perturbations in terms of dip using equations (2.10) and, finally, in terms 
of v using the definition (2.11). Note that, even though the relations (2.10) were derived by 
solving the constraints associated with the unperturbed action for the perturbations S2 So, 
corrections to these relations due to S2S1 do not contribute at first order in A. To show this, 
let us expand and B in powers of the coupling 



+ A0i + O(A 2 ), B = B + XB 1 + C(A 5 



(4.1) 



where 0o and Bo satisfy the unperturbed relations (2.10a) and (2.10b) respectively. By 
expanding the full quadratic action for the perturbations 82S [Sip, <fi, B] to first order in A we 
obtain 



S 2 S [Sip, 0, B] ^S 2 S [<V,0 O , B ] + A 



S(S 2 S 



k + x 



A=0 



S(S 2 Sq) 
SB 



Bx + 



A=0 



+ \S 2 S 1 [Sp,^B } + O{X 2 ). 



(4.2) 



However, the second and the third term vanish because they contain the unperturbed con- 
straint equations evaluated at 0o and Bq, which by assumption are solutions to the con- 
straints. Therefore, as long as we are interested in first order results, we can safely use the 
unperturbed solutions 0o an d -Bo given by equations (2.10). 

Expanding the leading action So to quadratic order in the perturbations, we obtain the 
free action (2.7). The additional quadratic terms stemming from Si must be appropriately 
contracted expressions containing partial derivatives of the perturbation variable v, the scale 
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factor a and the background field ipo. However, in the case of scalar perturbations, the field 
v can arise from fluctuations of the scalar field, 5ip = v/a, or from fluctuations of the metric, 



aV2e 



2v k 6Vo + ^ K - aHv k ) (5^5% + SVj) 



aV2e V» v 



(4.3) 



Thus, unlike the case of tensor perturbations, 8g^ u provides an additional tensor structure 
that can be used to contract derivatives. We now show that such contractions yield terms 
where the derivatives acting on v or a are contracted with r]^ v '. This means that the argument 
in the previous section can be applied to scalar perturbations as well, yielding essentially 
the same results. We note that terms which contain only fluctuations coming directly from 
the scalar field do not present this problem, and can be easily written as in equation (3.8). 

Let us first consider terms with only one factor of V^ u . In this case, V^ v can be contracted 
either with 77^, leading to V^r}^ = 2v, or with two derivatives <9 M ct, resulting in 



' dnd^a 
2 1 M 



dudd^a 



c^k^lV^ = 2{d^[a,<p ])v 



(4.4a) 

(4.4b) 
(4.4c) 



where, again, the square brackets in the last line mean that the derivatives can act on one or 
more factors of a or tpo. Thus, terms with only one factor of do not present any problem 
since, as anticipated, all the derivatives are contracted with the inverse of the Minkowski 
metric. 

Corrections which contain two factors of V^, and are not products of terms in (4.4b), 
can always be recast as 



2v 2 + 



k 2 



d„a 



d,,ad^a 



a a 2 



(4.5a) 



a 



d^ad^a 

-dyV a v 



d^a d IM ad ll a\ d u a 2 



a 

d u (v 



(4.5b) 



a \ a a* / a 

All the terms inside the square brackets become negligible in the sub-horizon limit, since 
their contribution is suppressed by an extra factor of 1/k 2 . The only term in which some 
derivatives are not contracted with the Minkowski inverse metric is the first one in equation 
(4.5b). However, the two derivatives with respect to conformal time result in a factor of k 2 
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which is precisely canceled by the extra factor 1/k 2 , and for all practical purposes such a 
term is equivalent to d^v d^v. 

Therefore, we have demonstrated that terms quadratic in the scalar fluctuations can be 
schematically written as in equation (3.8). The remainder of the analysis then proceeds 
as for tensor perturbations, and effective corrections to scalar perturbations are thus also 
sub dominant in the regime 

M 2 

H < M p and k ph < A ~ -f . (4.6) 

Before concluding, we would like to address again whether the operators that we have 
considered can be actually obtained from generally covariant terms. In the case of scalar 
perturbations, it is indeed possible to show — after further rather lengthy calculations — that 
the following family of covariant terms generates the kind of corrections shown in Table I, 

d-f = 0: R^(V^)V uV 

d-f = l: i^(V^)V,V 7 V> (4.7) 

d - f = 2 : (V a V p Rn (V Q V M <p) V P V U y? 

d-f = 3: (V Q V^) (V a V, <p) V^V, V 7 V 7 if 



In order to illustrate how this happens, let us consider for example the d — f = 1 term. It 
contains, among many other terms a factor 

a 2 BT d, 5<p d v 9 7 9 7 5<p d 4 d^a d u a d, v d v d^d~<v~ -% d"a d u a d,d v vUv + ... (4.8) 

where, in the last step, we have neglected a subdominant contribution in the short- 
wavelength limit. The last term in (4.8) indeed generates a correction proportional to 
H 2 k 2 h /Mp and it is suppressed by one factor of the slow- roll parameter. It is relatively easy 
to verify that the corrections generated by the other members of the family (4.7) have the 
same slow-roll suppression, which strongly supports the assumption we made in the context 
of tensor perturbations. 

V. CONCLUSIONS 

The connection, through cosmological inflation, between physics on the smallest scales, 
described by quantum field theory, and that on the largest scales in the universe is one of 
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the most profound aspects of modern cosmology. However, since inflation takes place at 
such early epochs, and magnifies fluctuations of such small wavelengths, it is important to 
establish the regime of validity of the usual formalism — that of semiclassical gravity, with 
quantum field theory assumed valid, and coupled to the minimal Einstein-Hilbert action — at 
those scales. 

On general grounds, we expect the canonical approach to break down at ultra-short dis- 
tances, where the operators that arise in an effective field theory treatment of the coupled 
metric- inflaton system become relevant. In this article we have calculated the impact of these 
higher-dimensional operators on the power spectrum at short wavelengths. In this way, we 
have been able to determine the regime in which the properties of the perturbations deviate 
from what is conventionally assumed. From a purely theoretical standpoint, these consid- 
erations are important if we are to understand the limits of applicability of cosmological 
perturbation theory. From an observational standpoint, cosmic microwave background mea- 
surements are becoming so precise that we may hope to use them to identify the signatures 
of new gravitational or field theoretic physics. 

Our analysis has focused on tree-level corrections to the spectrum. Because we have 
essentially considered all possible generally covariant terms in the effective action, we expect 
to have unveiled the form of all possible corrections that are compatible with the underlying 
symmetries of the theory. Loop diagrams may yield additional scale-dependent logarithmic 
corrections neglected in our analysis. In any case, our results indicate that cosmological 
perturbation theory does not apply all the way to infinitesimally small distances, k p h — > oo, 
and that, indeed, there is a physical spatial momentum A (or a physical length 1/A) beyond 
which cosmological perturbation theory breaks down. The scale at which perturbation 
theory breaks down at tree-level is 

M 2 

A~-^, (5.1) 

which, because of existing limits on the scalar to tensor power spectrum ratio [38] , is at least 
10 4 times the Planck scale. 

These results have significant implications for the impact of an effective cut-off A on the 
primordial spectrum of primordial perturbations, which typically is at most of order H/A 
[19]. Substituting the value of A we have found, we obtain corrections of the order of H 2 /M 2 , 
which are likely to remain unobservable [23] . This value of A also solves a problem that was 
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noticed in [39], namely, that in the presence of a Planckian cut-off, cosmological perturba- 
tions do not tend to decay into the Bunch-Davis vacuum (or similar states). In particular, 
to lowest order in perturbation theory, the transition probability from an excited state into 
the Bunch-Davis vacuum is significantly less than one for A = M p , but proportionally larger 
if A is given by (5.1). Ultimately, a large decay probability is what justifies the choice of the 
Bunch-Davies vacuum as the preferred initial state for the perturbations at scales below the 
cut-off, since, as we have found, our theories lose their validity at momentum scales above 
the spatial momentum A. 
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APPENDIX A 

In this appendix, we show how to integrate by parts every term of the form 

d 2n+m+p d 2 q +m+r y gto+p+r v ( A1 ) 

in order to express it as a linear combination of terms like 

d 2j+l [a,^]d 2m+l vUv (A2) 

plus, possibly, a term with no derivatives acting on v. Notice that, for notational convenience, 
we have defined □ = d^d^. Of course, if the index q (or s) in equation (Al) is not zero, 
we can easily integrate by parts 2q + m + r — 2 (2s + p + r — 2) times to get only terms of 
the form of that in equation (A2). Therefore, in what follows we only consider terms with 
q = s = 0. In this case, we can always integrate by parts an appropriate number of times to 
get only terms for which m = p. Thus, without loss of generality, we can restrict ourselves 
to considering terms of the form 

d 2n+ m+P ^ ^ grn+r y gp+r ^ (m = P ) . (A3) 
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The derivatives acting on v that are contracted with derivatives acting on a or (fo can be 
systematically eliminated by repeated integrations by parts: 

d 2n+m+P [a,Vo}d m+T vd P+r V d 2n+ra+(p-l) ^ Q rn+r y q (p-l)+r Q y (A4) 

+ i 9 2(n+i)+(m-i)+(p-i) [ fl)¥ , ]9(M+W o gM+W Vi 

where we have denoted equivalence up to integration by parts with ~. The first term on 
the right hand side can be cast in the form (A2) by integrating by parts (p — 1) + r times, 
while the second one is of the form (A3) with n and r (p and m) increased (decreased) by 
one. By iterating this procedure, we eventually obtain terms of the form 

d 2n [a, Vo ]d r vd r v , (A5) 

where now n and r have changed. Again, we can integrate by parts and obtain 

d 2n [a, ifo] d r vd r v~-d 2n [a, <p ] d r ~ l vd r ~ l □ v + ]-d 2{n+1) [a, <p ] 9 (r " 1} v d (r " 1} v . (A6) 

The first term on the right hand side can be re-written as (A2) after r — 1 integrations 
by parts, while the second term has the form (A5) with n (r) increased (decreased) by 
one. Thus, by repeating this procedure we obtain many terms of the form (A2) and we are 
eventually left with a term without derivatives acting on v. This completes our proof. 
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